rjc mixing effects on charmonium and B meson decays 
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We include the rjc meson into the rj-ri'-G mixing formalism constructed in our previous work, where 
G represents the pseudoscalar gluball. The mixing angles in this tetramixing matrix are constrained 
by theoretical and experimental implications from relevant hadronic processes. Especially, the angle 
between rjc and G is found to be about 11° from the measured decay widths of the rjc meson. The 
pseudoscalar glueball mass ma, the pseudoscalar densities mqq^s3,cc and the U{1) anomaly matrix 
elements associated with the mixed states are solved from the anomalous Ward identities. The 
solution ma ~ 1.4 GeV obtained from the rj-ri'-G mixing is confirmed, while m^g grows to above the 

■ pion mass, and thus increases perturbative QCD predictions for the branching ratios Br{B — > rj'K). 
We then analyze the T^c-mixing effects on charmonium magnetic dipole transitions, and on the 
B — >■ rf'^Ks branching ratios and CP asymmetries, which further improve the consistency between 

■ theoretical predictions and data. A predominant observation is that the ?7c mixing enhances the 
CO ■ perturbative QCD predictions for Br{B — > rj' K) by 18%, but does not alter those for Br{B — s- rjK). 

The puzzle due to the large Br{B — >■ rj'K) data is then resolved. 

Oi 
D 

I. INTRODUCTION 

^ ' It has been known that the gluonic and charm contents of the light pseudoscalar mesons r] and rj' may have a crucial 

■ impact on studies of many hadronic processes, such as 77^'-* electromagnetic (EM) transition form factors, 77*^'^ radiative 
] decays, and charmonium and heavy-flavor decays into 77*^'^ . For a recent review on the gluonic effects in EM transitions 

■ and in weak decays of charm and beauty hadrons, see [ij. It remains a puzzle that theoretical predictions for the 
B — > r]'K branching ratios are usually lower than data, even after taking into account the 77-77' mixing 0,3- Hence, 
it has been conjectured that the gluonic content of the 77' meson plays a role in accommodating the large branching 

A gluonic content of the 77^'^ meson has been inferred from data of the radiative decays P — )• and 
[To| and from the charmonium decays J/^ ^ VP [HI, However, the gluonic contribution to the 
B — )• 7;^'^ transition form factors was parametrized and tuned to fit data in the QCD factorization approach [Tsj and 
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in the soft-coUinear effective theory [1J|, so no conclusion on its importance could be drawn. This contribution was 
calculated explicitly in the perturbative QCD (PQCD) approach [T^l using the gluonic distribution amplitudes of the 
' 77*-'^ meson from [16| and in QCD sum rules [17], and it was found to be small. 
5^ , The charm content of the 77^'^ meson has been introduced through the ri-ri'-rjc mixing [Tsl. figj. This formalism was 
extended to the tetramixing among the tt, ry, 77', and rjc mesons recently 20], whose parameter set, including the 
mixing angles and the hadronic parameters in the light-front constituent quark model, was determined by a fit to 
data of relevant meson transition form factors. An intermediate question is whether this charm content affects the 
B rj^'^K branching ratios [2l[ and their CP asymmetries [l^, whose measurement might reveal new physics signals. 
A potential deviation has been detected between the mixing-induced CP asymmetries in the tree-dominated decays 
B — > J /ipKs and in the penguin-dominated decays B — >■ 77'^^^ . Whether this deviation can be interpreted as a signal 
of new physics depends on how large the tree pollution in the latter is. Though the ri^'^-r]c mixing is small (23l. |24| . 
there is lack of quantitative estimate of its effect. It is then worthwhile to examine whether the large B — rjcK 
amplitudes are able to compensate for the tiny mixing, and to give a sizable impact on the B — r]^'^K decays. 
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The above subjects demand complete and precise understanding of the gluonic and charm contents in the r]^'^ meson. 
We have set up the rj-rj'-G mixing in our previous work [25|, where G denotes the physical pseudoscalar glueball. This 
mixing was implemented into the equations of motion from the anomalous Ward identities, that connect the vacuum to 
77, rj' , and G transition matrix elements of the divergence of axial- vector currents to the pseudoscalar densities and the 
C/(l) anomaly. Solving these equations, the pseudoscalar glueball mass niG was expressed in terms of phenomenological 
quantities such as the 77, 77' masses, the decay constants, and the mixing angles. With the mixing angles measured 
from the (f> — 77/, 77;' decays by KLOE mc w 1.4 GeV has been deduced, suggesting that the 7;(1405) meson is an 
interesting pseudoscalar glueball candidate [25-27]. However, the solution for the pseudoscalar density mqq associated 
with the rjq component of the rj^'^ meson is still as low as its conventional value 777^^ « m^r ~ 0.14 GeV. It has been 
postulated that a larger mqq « 0.2 GeV can enhance the B ^ rj' transition form factors, and thus the B rj'K 
branching ratios significantly [28j . Following this vein, it was pointed out that the introduction of decay constants, 
suppressed by the Okubo-Zweig-Iizuka (OZI) rule |29l - l3ll |. into the equations from the anomalous Ward identities can 
increase mqq [isj . 

We are motivated to formulate the tetramixing among rj, rj' , G, and 77c, and to investigate its impact on charmonium 
and B meson decays in this paper. The mixing with the pion is not considered here under the isospin symmetry. As 
a consequence of the mixing, the 77c meson contains a gluonic content, that modifies the QCD calculation of its decay 
width. The fit to the observed 77c decay width determines the additional mixing angle (pg ^ 11° between G and 77c. 
We shall explain that the gluonic content of the ?7c meson further improves the calculations for the decay widths of 
the charmonium magnetic dipole transitions J/V', ip' — 777c in association with the unquenched mechanism proposed 
in [3^ Issf. and renders theoretical predictions in better agreement with data. Together with the other mixing angles 
fixed in l25ll , we construct the tetramixing matrix, which implies the charm content of the rj' meson consistent with 
that in [20|. Implementing the tetramixing into the equations from the anomalous Ward identities, we solve for the 
pseudoscalar glueball mass, the pseudoscalar densities, and the U{1) anomalies. It is found that the inclusion of the 
77c mixing does not alter our prediction for the pseudoscalar glueball mass, but increases the pseudoscalar density mqq 
to above the pion mass, even in the absence of the OZI-suppressed decay constants. 

Moreover, the charm content of the 77''^ meson allows the B — >■ rj'^'^K decays via the B — rjcK channel. Simply 
adopting the amplitudes evaluated in the PQCD approach at the next-to- leading-order (NLO) accuracy [ll, HI], we 
estimate the 77c mixing effects on the B 'ii'--'^K decays. It will be demonstrated that the additional tree contribution 
from B — rjcK increases the B — rj' K branching ratios by about 18%, but does not change the B — 77^^ branching 
ratios. Combining the mechanisms from the larger mqq and the charm content, the puzzle due to the large B — >■ -q' K 
branching ratios is resolved. On the other hand, the charm content of the t/''' meson has minor infiuences on the 
direct and mixing- induced CP asymmetries in the B — rj'^'^Ks decays. Nevertheless, we do see the modification 
toward accommodating the deviation between the measured mixing- induced CP asymmetries in the _B — J/ ij^Kg and 
B — > r]' Ks decays. The 77c mixing effects on the CP asymmetries in the B — > rjKs decays are also negligible. 

We set up the r]-!]' -G-rjc tetramixing formalism in Sec. II, which contains one more mixing angle 0q between G 
and 77c compared to the rj-rj'-G mixing. The angle 4>q is then determined, and the resultant gluonic and charm 
contents of the 77*^'^ meson are compared with those obtained in the literature. In Sec. Ill we solve for the pseudoscalar 
glueball mass, the pseudoscalar densities, and the U{1) anomalies appearing in the mixing formalism, and discuss 
their phenomenological implications. The ?7c mixing effects on charmonium magnetic dipole transitions and on the 
B — > rj^'^K decays are also investigated. Section IV contains the summary and comments on other works, that present 
observations different from ours. 

II. n-v'-G-Vc MIXING 

In this section we formulate the rj-rj'-G-rjc tetramixing, determine the involved mixing angles, and implement the 
mixing into the equations of motion from the anomalous Ward identities. 

A. Mixing matrix 

We combine the Feldmann-KroU-Stech (FKS) formalism for the rj-rj'-rjc mixing [l^, [l^ and for the rj-rf-G mixing 
[25} , in which the conventional singlet-octet basis and the quark-flavor basis qq = {uu + dd) / \/2 and ss |3^] , labeled 
by rjq and 77^, respectively, were adopted. We further introduce the unmixed glueball state g and the unmixed heavy- 
quark state 77Q. Let the matrix U^i (C^i4, U12) represent a rotation with the 3-4 (1-4, 1-2) plane flxed. It is natural 
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to first mix those singlet states rji, g, and rjQ, and then mix 771 and rjs to form the physical states 77, 77': 

= C/34(0)t/i4(0G)t/l2(0Q) 




\9) 

Wvq) 



(1) 



with the rotational matrices 



/cose - sine 

sin6' cos6' 
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V 1 





1 











COS0G sin(/>G 
— sin (f)G cos (f)G 
1 



(2) 



That is, we assume that the octet state rjs does not mix with the glueball, and that the heavy-flavor state mixes with 
the pseudoscalar glueball more dominantly than with the 771 state. 
The octet and singlet states are related to the flavor states via 



/ \Vs) 



\9) 

Wvq) 



where 0^ is the ideal mixing angle with cos6'i — yl/S and sin^?i = y^/S, i.e. 
transformed into the physical states through the mixing matrix 




(3) 



— 54.7°. The flavor states are then 



U{e,cbG,cl>Q) = U34{e)Uli{cj,G)Ul2{^Q)U34{e^), 

/ c^c^i — seapGsei —cesei — s 



secei + ceapGsei 
s(t)Gse.i 




s^s^i + cec(f)Gcei 
-scjiGcei 





~ses(j)GC<j)Q 
ces(f)GC(t)Q 

C(t)Gc4>Q 

-s(t>o 




(4) 



limit, 



with the notations c0 = cose and s0 = sin0. Equation ([4]) approaches the mixing matrix in [25| in the 
and the 77-77' mixing matrix fl§L [l9| in the (f>Q, (pG — > limit. 

As stated in the Introduction, we have assumed isospin symmetry, i.e. no mixing with 7r°, and ignored other possible 
admixtures from radial excitations. The decay constants /g, fs and fc are defined via the matrix elements 



)|s7'^75s|77.(P)> 
|c7^75chQ(P)> 



-ifcP^, 



(5) 



for the light quark q = u or d. The other decay constants of the rjq, rj^, and 77o mesons and of the unmixed glueball, 
which are suppressed by the OZI rule, can be introduced in a similar way [l5|: 



mi^l,q\Vs{P).9{P),VQ{P)) = -^fl,,,P^, 

{0\rsr-/5sMP),9{P),VQiP)) = -^Ila,cP^. 
mri,c\ij,{P),ij,{P),g{P)) = -ifls.gP^. 

The decay constants associated with the 77, 77', G, and 77c physical states in the following matrix elements 

mri5q\ri{PU{P),G{P).ric{P)) = g,,. 



(6) 



(0|s7^75sh(^),^'(^),G(P),77,(P)) - -7/,^, 
mriMP).l'{P)MP).Vc{P)) - -ifl^> 



G.Vc 



(7) 



4 



are related to those associated with the ijq, r]s, g, rjQ states through 



/ fv fv fv \ 

fs 

fk fs fc 

JG JG JG 

\ f f f J 



( h n \ 

J s J s J s 

f? fs fc 

Jg Jg Jg 

\Pc Ic fc) 



(8) 



We sandwich the equations of motion for the anomalous Ward identities 

a. — 



An 



(9) 



between vacuum and \ri), \ri'), |G), and {rjc), where 



are the quark masses, Gfj,i, is the field-strength tensor and 



G^"^ is the dual field-strength tensor. Following the procedure in [15|, we derive 



M, 



qsgc 



^u\e,4>G,^Q)M^u{e,4>G,^Q)l 



in which the matrices are written as 



M. 



qsgc 



/ mlq + V2Gg/.f, m% + Gg/.fs m% + Gg/U\ 
m% + ^Gs/f, ml + GJfs m^^ + Gs/fc 



m 



qg 







( < , 








^Ggl]g m%+Gglfs 711^ + G g / 



V2GJf, 






to2. 



V 



\ 





m^J 



Gc/fs ml, + GJfJ 

( 1 PqlJs f^/fc\ 



J = 



n/fq 



nifc 



fg / fq fg/ fs fg I fc 
\fllfq fllfs 1 / 



(10) 



(11) 



(12) 



with the ?7, 77', G, t]c meson masses ?ti^_^/_g_^^, and the abbreviations for the pseudoscalar densities and the [/(I) 
anomaly matrix elements 



qq.qs.qg.qc 
•,2 

"sq.ss^sg^sc 
cq.cs,cg,cc 



V2 

fq 



(0|m„m75U + mddij5d\r]q, r]s,g, r]Q) 

0\msSi^5s\r]q,r]s,g,r]Q), 

Q\mcCij5c\'qq,'i]s,g,'nQ), 



G 



q,s,g,c 



(0|a,GG/(47r)|?7,,77„g,77Q). 



(13) 



In the following analysis we neglect all the OZI-suppressed decay constants defined in Eq. (|6|) according to the FKS 
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scheme. Equation p2)) then leads to the following equations explicitly 

K%c = -ml{c9c0,-s9c(^GsO,){c9s9, + s9ccj)Gce,) 

+mfy{s9c9.i + c9c(j)Gs9i){-s9s9^ + c9c(j)Gc9i) + m^(j{s(t>Gfc9is9i, 

Ml!gc = -m'^^{c9c9,- s9c<i)Gs9i)s9s(j)GS(t)Q 

+m'^,{s0c9i + c9c(l)Gs9i)c9s4)GS(l)Q - m%c(l)GS(t)Gs9iS(l)Q, 

qsgc qsgc^ 

Mqsgc - rnlic9s9, + s9c(j)Gc9,f + ml,{~s9s9, + c9c(j)Gc9,f + ml{s(j)Gc9,f , 
Mqsgc = m^{c9s9t + s9c(t)Gc9i)s9s4>GS(t)Q 

+m'^,{-s9s9t + c9c4>Gc9i)c9s4>GS(t)Q - m%c(l)GS<j)Gc9is4>Q, 
^ifsgc = -ml{c9c9,- s9c4>Gs9,)s9s4>Gc4>Q 

+mfy{s0c9i + c9c(t)Gs9i)c9s4)Gc4>Q - mQC(j)GS(j)Gsdic4>Q, 
Mqsgc = m^{c9s9i + s9c(l)Gc9i)s9s4)Gc4)Q 

+m'^,{-s9s9t + c9c4)Gc9i)c9s4)GC(t>Q - ■m%c4>GS(t>Gc9ic4>Q, 

^"^qlgc = m'^{s9s(f>Gfc(f>QS(f>Q +m'^,{c9s(f>Gfc(t)QS(f>Q 

+m%{c(t)Gfc(j)QS(t)Q - m^j^c(t)QS(f>Q, 
M^lg^ = -m'^{c9c9i - s9c4)Gs9i)s9s4iGS(l)Q 

+m'^,{s9c9i + c9c4'Gs9i)c9s4>GS(t>Q - m'^Qc4>Gs4>Gs9is4>Q, 
Mfsgc = ml{c9s9, + s9c(t)Gc9,)s9s(t>GS(t)Q 

+m'^,{~s9s9i + c9c4>Gc9i)c9s4iGS(t)Q - m%c(t)GS(t)Gc9is4}Q, 
Mq!gc = ml{s9s^GS(l)Qf + ml,{c9sc^GSC^Qf + ml{cc^GS<t>Qf + ml^{c^Qf . 

B. Mixing angles 

There is already extensive discussion on the determination of the mixing angle 9 in the literature, whose value still 
varies in a finite range. For example, —17° < 9 < —11° has been extracted in [13], assuming the presence of the 
gluonic content in the r]' meson. We choose 9 = —11°, which corresponds to a sizable gluonic content [13], and is 
appropriate for our choice of below. The value of 4>g varies in a wide range, depending on the parametrization of 
the mixing matrix, experimental inputs, and fitting procedures [25j . Even its central value can take a number between 
10° < ^ 30°, such as (t>G = (12 ± 13)° in To] and (t)G = (33 ± 13)° in We take the value (/>g = 12°, which was 
also considered in [25| . 

The last angle (j)Q can be determined by the ijc total width and the rjc — > 77 decay width. As a cc(^So) state, the 
value of Ftot is quite large among the charmonia below the DD threshold. During the past few years, the experimental 
results for th e ?7p total width and the two-photon decay branching ratios vary drastically. The Particle Data Group 
(PDG) 2008 [If listed Ttot = 26.7 ± 3.0 MeV and Br{T]c 77) = (2.4^ J;^) x lO"'^, respectively. In contrast, the 
PDG2010 ^ presents Ftot 28.6 ± 2.2 MeV and Br{T]c -> 77) (6.3 ± 2.9) x 10"^. The BESIII Collaboration 
measured the ijc total width recently in ip' — ^ jrjc and found Ftot = 32.0±1.2±1.0 MeV [i^. The interesting tendency 
is that the total width of ijc becomes broader than the previous measurements. Such a change favors the scenario that 
the rjc has a glueball content in association with the cc(^So) component as analyzed below. One also notices the much 
smaller Br{ric — 77) listed in the PDG2010 than in PDG2008, which may be caused by the old total width data 
employed in the extraction of the two-photon branching ratio. Since the averaged partial decay widths for rjc — > 77 
are unchanged in PDG2008 and PDG2010, i.e. T{r]c — >■ 77) = 6.7to;8 keV, it might be more appropriate to adopt the 
experimental data for the two-photon partial width instead of for the branching ratio in the following analysis. 

The strong decay of cc(^So) annihilation via two-gluon radiation can be related to its two-photon decay. To lowest 
order, one has 

r(^^o ^ 99) J^^9/a,y 

F(i5o^77) " 9e4a2 8 \ae J ' ^ ' 



(14) 

(15) 

(16) 
(17) 
(18) 

(19) 

(20) 

(21) 

(22) 

(23) 

(24) 
(25) 
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where Be = 2/3 is the charge of the c quark, and Us and ae are the strong and EM couphng constants, respectively. 
With T{^Sq ^ gg) ~ Ttot = 28.6 ± 2.2 MeV from PDG2010 [H and T{^So 77) = 6.7^°;^ keV, we extract 
0.41 - 0.49 {as ~ 0.72 - 1.2 if one adopts the branching ratio Br{'qc lll^ (6.3 ±2.9) x 10"^ from PDG2010 
[39|), which is much larger than the running coupling as{m,c) — 0.24 ^ 0.26 [39|. Even with the first-order QCD 
correction [4l[ taken into account: 

r(i5o->.g5)^9^a,V(l + 4.8H.) ^^^^ 



r(i5o^77) 8 \aej (1-3.4^)' 

we still have larger values — 0.28 ~ 0.33. 

As defined in Eq. (jl]), the rjc meson has the wave function 

\Vc) =-s(l)Q\g) + c^Q\ilQ)- (28) 
The strong decay amplitude can be parametrized as 

{99\V\^c) = ~s^Q{9g\V\g)+c,pQ{gg\V\j^Q) 

^^ + ccj^Q]{99\V\m). (29) 



where V denotes the potential for annihilating the cc states and the glueball into two gluons. Since the glueball 
does not pay a price for coupling to the \gg) state, we have the ratio {gg\V\g) / {gg\V\'qQ) — i^jas given by the 
empirical gluon power counting [42l|. Therefore, the mixing provides a correction to the strong decays different from 
the correction to the EM coupling — (c(/)Q)^ag. Consequently, one has, to the leading QCD correction, 



r(^c ^ gg) _ 9 {-^^Q + (1 + 4.8^) 
r(77e^77) 8 al (1-3.4^)' 



(30) 



with t(^Q = tan^Q. Requiring ctg — as{m^ ~ 0.25, we deduce ~ —1.4° or 10.4° (0q ~ —6.9° or 15.6°, if 
Br(?7c -> 77) = (6.3 ± 2.9) x 10"^ is adopted). If inserting the recent BESIII value Ttot = 32.0 ± 1.2 ± 1.0 MeV [13], 
we extract 0q ~ —1.7° or 10.8°. It is noted in advance that the data of the radiative decays J/V', V'' ~^ IVc favor the 
positive solution as illustrated in Sec. Ill C. Based on the above determinations of (j)Q, we adopt (/jq = 11° in this 
work for the purpose of estimation. 

We should mention a recent calculation of the ^Sq decays into light hadrons and two photons in the framework of 
nonrelativistic QCD (NRQCD) [43], in which the 0{asV^) corrections were found to be crucial for accommodating 
the observed rjc — >■ gg and 77 widths. The involved long-distance matrix elements for the rjc meson still suffer large 
uncertainties, and their values were obtained by data fitting in (43| . Hence, it is a fair comment that the present data 
cannot distinguish the NRQCD results with the 0{asV^) corrections from the mixing scenario with a small glueball 
component in the rjc meson. Meanwhile, it should be realized that the presence of the glueball component would 
affect the total and EM decay widths of the rjc meson differently, and thus make an impact on their ratio at leading 
order. 

With the angles 9 = —11°, (f>G = 12° and c/)q = 11°, the mixing matrix in Eq. (j4]) is explicitly written as 



U 



I 0.720 -0.693 0.039 0.008 

0.673 0.711 0.200 0.039 

-0.170 -0.120 0.960 0.186 

V -0.191 0.982, 



(31) 



Compared to the parametrization for the ry-ry'-ryc mixing [3, [l9|, the role of the small angle — —1° ±0.1° has been 
played by our (t>G4>Q- The matrix for the rj-ri'-rjc mixing is given by [l8| 

\ri) = 0.77]77,)-0.63]r;,)- 0.006]ryQ), 
]??') = 0.63]77,)+ 0.77]r;,)-0.016]r;Q), 

]77,) = 0.015]ry,)+0.008]r;,) + ]r,Q), (32) 

in which the charm content in the rj^'^ meson has a sign opposite to that in Eq. (1311) . The fit to the data of 
r{J/il' — ^T]') /T{J /ijj — 777c) [3 actually cannot discriminate the sign of these coefficients. Another difference 
appears in the rjq^s components of the rjc meson. Since these are small components, there is no inconsistency between 
the results in Eqs. ^ and 
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The charm contents of the rj and 77' mesons have the same sign in the FKS scheme. This feature differs from the 
parametrization for the TT-rj-rj'-ijc tetramixing based on the group decomposition 50(4) — SO{3) (g) 5*0(3) ;20i] . 

0.9895 0.0552 -0.1119 0.0342 

I -0.1082 0.8175 -0.5614 -0.0259 . , . 

Unrir,'^, - | 0.0590 0.5696 0.8160 0.0452 ' ' 

-0.0395 -0.0065 -0.0478 0.9960 

where the charm contents of the rj and 77' mesons are opposite in sign. A careful look reveals that these matrix 
elements are small due to the destruction of large numbers, so they are sensitive to experimental inputs. Varying 
the inputs slightly, one could get the matrix elements of the same sign in The rj meson has a pion component 
—0.11, while the rj' meson has a smaller pion component 0.06 as indicated in Eq. (1331) . Therefore, it is appropriate to 
compare the charm contents of the rj' meson in Eqs. (|5T|) and both of which are about 0.04. 

III. 7?c-MIXING EFFECTS 

In this section we solve Eq. (|10l) from the anomalous Ward identities, and then investigate the phenomenological 
impacts from the ?7c mixing on charmonium magnetic dipole transitions and on the B — ^ rj^'^K decays. 

A. Decay constants 

To obtain the decay constants defined in Eq. (|5]), we need the inputs of fq, fs and fc for the fiavor eigenstates. 
The value of fq is close to the pion decay constant Z^, such as fq — (1.07 ± 0.02)/7r extracted in [l^. Our analysis 
indicates that the variation of fq has almost no influence, so we simply set it to fq = f-K = 131 MeV. The value of fs is 
more uncertain, for which the result fs/ fq ~ 1-2 ~ 1.3 also from [18| is employed. Since the ratio fs/fq gives a more 
significant effect, we shall examine how our outcomes depend on its variation. For fc, we adopt the value fc = 487.4 
MeV 0]. Ignoring the terms suppressed by the OZI rule, Eq. ^ leads to 



p 
J77 


fv] 












ffr 









fs 





fc 


fc 











fl 


ft) 




V 





fc 



113 -90.8 3.69 

106 93.1 19.0 

-26.7 -15.7 91.0 

478 



(34) 



The larger decay constants are close to those appearing in the literature [18], and in the perturbative calculations 

of the B — ^ rj^'^K decays [H, ^3, H^l- It is natural that f^^ is almost the same as fc- Those decay constants /g'*''^ 
may have phenomenological applications, after the pseudoscalar glueball is identified. Because of = —11°, we have 
the ratio f^/ f^, = - tan 61 = 0.24, which is similar to -0.006/(-0.016) = 0.38 in [11 El, but has an opposite sign 
to —0.0259/0.0452 = —0.57 in [2^. Note that the sign of ffy is still not certain, which was found to be positive in 
[21I [4^ , but negative in [11 [11 . According to [1^ , a positive /^j,, would enhance the 77*^'-' transition form factor in 
the light-cone PQCD approach. Nevertheless, the enhancement will be compensated by varying other parameters, 
such as the constituent quark masses, so that the corresponding data can still be accommodated in their theoretical 
framework [46| . Our result /^/ = 19.0 MeV, arising from the phenomenological determination, is approximately equal 
to 21.9 MeV in [2^, and lower than 50-180 MeV computed from QCD low energy theorem [21]. The similarity to 
the value in [20j is nontrivial, since we have constructed the mixing matrix via the angles Q, (po and (j)Q, which were 
determined in a different way. It has been conjectured [l^l that the value of /,'^, was overestimated in [2l|. Hence, 
we disagree with their speculation that the charm content of the 77' meson alone can exhaust the large B — > rj'K 
branching ratios. 

Our result f^, = 19.0 MeV is larger than —(6.3 ± 0.6) MeV in magnitude in [l^, and about 8 times larger 
than 2.4 MeV extracted from the data of Br{J/ilj — )■ jr]')/Br{J/tp — 777c) [IJ. As stated above, the data of 
Br{J/ifi — > ■f?]')/Br{J/^ — > jrjc) cannot fix the sign of f^y in [ll|. Besides, the smaller value in [ll| is attributed to 
the tiny 77^ mixing, which is a consequence of the anomalous Ward identities in Eq. ^ without the OZI-suppressed 
terms. Including the OZI-suppressed terms, a larger range for is expected. The analysis in [H, [13], together 
with that in [48l . l49j which also concluded a small 7;'-77c mixing, were performed in the rj-rj'-i^c mixing formalism. 
That is, only a single channel J/tp — ^rjQ — )■ 777' contributes to Br{J/il) — 77;'), which is perhaps too strong of 
an assumption. With the more general rj-rj'-G-rjc tetramixing, an additional channel J/ip — > 71? — > 777' exists. For 
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the purpose of ilhistration, we regard 77(1405) as the pseudoscalar glueball ^2^, and assign a finite fraction of the 
J/V' 777(1405/1475) amplitude to J/t/" 777(1405). Employing the data in PDG2010 SQi], it is easy to find that the 
above two channels, multiplied by the mixing matrix elements 0.039 and 0.2 in Eq. pi|) . respectively, are comparable 
to each other. If a destructive interference occurs between the two channels, a larger 77^-77^ mixing will be allowed. 
On the other hand, some theoretical estimates on also implied a small 77'-77c mixing [Tsl. [23l 47 [, in which the 
annihilation of the cc pair into two gluons, followed by their combination into the 77' meson, was considered. Strictly 
speaking, what they investigated is the extrinsic charm content of the 77' meson, while the 77'-77c mixing results in the 
intrinsic charm content [20]. Therefore, their observation does not contradict to ours, and supports that the intrinsic 
charm content analyzed in this work might be more relevant. 



B. Glueball mass, Pseudoscalar densities, and U{1) anomalies 

Given the decay constants fg, fs and /c, we solve Eq. (|10p in order to get an idea of the magnitude of the pseudoscalar 
gluball mass, the pseudoscalar densities, and the U{1) anomaly matrix elements. There are 12 equations from the 
4x3 matrix equation in Eq. ([TU)) . Even having dropped all the OZI-suppressed decay constants, there are still too 
many unknowns listed above. Hence, we intend to drop the OZI-suppressed pseudoscalar densities, except 777^^, 
and rn^g, whose values can serve as a check of their scaling behavior in the large Nc limit [2^, Nc being the number of 
colors. Then one of Eqs. (^01), (ED), (IMl), and (Ull) becomes redundant, since both the ratio of Eq. over Eq. (PTI) . 
and the ratio of Eq. ([^ over Eq. ([M]) give 

-rn^icOcOi - s9c<j)Gs0i)s6 + m^, {sOcOi + c9c4)Gs9i)cO - mQC(t>Gs9i \f2js 

m'^{cdsdi + s0c(l)Gc6i)se + m'^,{-s9s9i + c9c(j)Gc9i)cd - m'^c(l)GC0i ~ fq ' ^^^^ 

which is identical to the formula derived in (25j . Eventually, we solve for 11 unknowns, which include the pseudoscalar 
glueball mass uig, 6 pseudoscalar densities rn^^, ™L, "^cg, "^cg, and 77z^^, and 4 [/(I) anomaly matrix elements 
Gq, Gs, Gg, and Gc- 

According to Eq. (fT4)) . 771^^ is expressed as the difference between the right-hand side of Eq. (fT4|) and the anomaly 
matrix element Gg/ fq. Equation (|15|) implies that Gq is proportional to /s, so the above difference strongly depends 
on the ratio fs/ fa [15|. Equation psp shows that the inclusion of the 77c mixing does not affect much the solution 
niG ~ 1.4 GeV [25| . The value of mc is insensitive to the uncertain mixing angle 0g also, because of c0g ~ 1 for a 
small (j)G- Therefore, ruG is only sensitive to the ratio fs/ fq, and we shall consider the variation of this ratio below. 
The solutions corresponding to fs/ fq = 1.2 and 1.3 are collected in TableHl It is seen that only mc and 777^^ depend on 
the ratio fs/ fq, and other quantities are relatively stable. The results of niG, ^tiLj Gq, Gs, and Gg are similar to those 
from the rj-rj'-G mixing formalism [25|. Namely, 77i^„ in Table H] respects the leading Nc relation m^^ = 2m^ — trJ, and 
we do not observe the enhancement claimed in [50|. Note that a larger pseudoscalar glueball mass (> 2 GeV) has been 
postulated in a dynamical analysis of the mixing in the pseudoscalar channel [sTf. However, if their assumption on the 
meson couplings is relaxed, a lower mass can be attained. Our solutions respect the hierarchy [tti^^I 3> \'m'^g\ ^ Ito^^^^sI 
in the large Nc limit. The values of tti^^ are consistent with the relation 777^^ w 777^^ [l^, but the magnitude of Gc is 
a bit smaller than that in [l8| . 

It has been shown that the B ^ rj' transition form factors and the B — >■ -q' K branching ratios are sensitive to 
the pseudoscalar density 777^^ psj , which defines the normalization of the two-parton twist-3 distribution amplitudes 
for the rjq state. Its value is usually assumed to be the pion mass, mqq ~ mrr, with which the branching ratios 
Br(5± -J> v'K^) « 51 X 10"^ and Br{B° ■q'K") sa 50 x 10"^ have been obtained in NLO PQCD [34]. These 
results are stih lower than the data Br{B^ -q' K^) « (71.1 ±2.6) x 10"^ and Br{B° -q' K°) « (66.1 ±3.1) x 10"^ 
[5^ . It was then demonstrated that mqq can be increased up to 0.2 GeV by introducing the OZI-suppressed decay 
constants /| and /I into the equation for the q-q' mixing [l^. With the enhanced niqq, it is likely to explain the large 
Br{B — >■ r}' K) data with larger B ^ q' transition form factors. Since and /| are free parameters, whether mqq 
can reach 0.2 GeV is not conclusive. Table U indicates that the 77^ mixing, which receives a phenomenological support 
as elucidated in Sec. II B, can enlarge m^^ for fs/ fq = 1-3 by a factor 2, from 777^^ ~ 0.012 GeV^ in Eq. (29) of 
to 0.023 GeV^. Note that the values of the mixing angle 9 for Eq. (29) in (25| and for Table |T] are different. However, 
777g„ does not much depend on the variation of 9. 



C. Charmonium Ml transitions 



There has been a long-standing puzzle from the magnetic dipole (Ml) transition J/ip — >■ 777c. In contrast to the 
success of the nonrelativistic potential models in offering an overall good description of the charmonium spectrum. 
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TABLE I: Solutions corresponding to fs/fq ~ 1-2 and fs/ fq ~ 1-3. 





ma (GeV) 


2 2 
^qq ^^ss 






2 2 


(GeV^) 


Gq 


Gs 


Gg 


Gc (GeV^*) 


f. 1.2 


1.519 


0.067 0.443 


-0.156 


-0.092 


-1.197 5 


.648 


0.053 


0.031 


-0.023 


-0.004 


U 1.3 


1.376 


0.023 0.457 


-0.149 


-0.081 


-1.283 5 


.631 


0.056 


0.030 


-0.016 


-0.003 



the predicted Ml transitions between the vector and pseudoscalar charmonia appear to have significant discrepancies. 
Namely, the predicted partial d ecay width T^^{J/%p — ^ 777c) ~ 2.4 ^ 2.9 keV [H, [sj] is obviously larger than the 
experimental data in PDG2010 [39], Br{J/ij 777^) = (1.7 ± 0.4)%, i.e., T{J/^ -frjc) = (1.58 ± 0.37) keV. The 
PDG2010 value is mainly weighted by the CLEO data Br{J/ij ^ 777^) = (1.98 ± 0.09 ± 0.30)% [H. Although 
the CLEO data bring the experimental and theoretical values closer, the discrepancy remains nontrivial after taking 
into account another channel tp' •yrjc. the potential models predicted r^^{tfj' — )■ 77/c) — 9.6 ~ 9.7 keV, i.e., 
Brill}' -J> 777c) ^ 3.2%, while the CLEO measurement gives Br{ip' -> 777^) = (4.32 ± 0.16 ± 0.60) x lO'^ 

Theoretical efforts of stu dying t he charmonium EM Ml transitions in the framework of nonrelativistic multipole 
expansions can be found in [56l463l |. Recently, a nonrelativistic effective field theory was applied to J/ip — >■ 777c [62l|, 
in which the radiative decay width T{J/il) 777c) = (1.5 ± 1.0) keV up to correction of 0{v1/rn^) was obtained with 
a rather large uncertainty. This approach becomes unreliable in ip' — 77yc, because the cc pair cannot be treated as 
a weakly bound system anymore. The lattice QCD calculations of these processes were reported in [g^, H^. In the 
quenched approximation the result for J /ip — )■ 777c turns out to be in agreement with the potential models, while that 
for ip' — 777c is much smaller and compatible with the data within uncertainties. That is, the quenched lattice QCD 
does not resolve the puzzle completely. 

A possible resolution arises from the accommodation of open threshold effects as an unquenched mechanism in the 
charmonium Ml transitions [s^ [s^. Because of the presence of the open DD threshold, the transition ip' — > 777c 
would experience more influence from the DD threshold, which naturally lowers the partial decay width predicted 
by the potential models. In contrast, the open threshold effects on J/tA — )■ 77;c are relatively small, since the mass of 
J/ip is located rather far away from the DD threshold. However, the uncertainties from the unquenched effects are 
significant as shown in _32, 33], so there is still room for the glueball- r/c mixing mechanism in the present experimental 
and theoretical situations. 

Starting with Eq. (pS)) . we express the quenched Ml transition amplitude as 

{Vc\Heni\J/llj,1p') = -S(j)Q{g\Hera\J/i^,i^') + C(t)Q{r]Q\Hera\J/llj,1p') 

- {~S(j)Q— + C(t)Q){l]Q\Hem\J/^J,'lp') , (36) 

where {rjQ\Hp„i\J is equivalent to the potential-model Ml transition amplitude for the J/ip^ ijj' mesons, and 
the gluon counting rule has been implemented in the second line. One immediately notices that in order to lower 
the Ml transition partial width, a positive (j)Q is favored. Given (j)Q ~ 11°, the quenched Ml transition partial 
widths are lowered by about 7%. That is, the glueball-r/c mixing does improve the overall consistency between the 
potential-model predictions and the data for the charmonium Ml transitions. 



D. v'^'^Ks decays 

As stated in the Introduction, a potential deviation has been detected between the mixing-induced CP asymmetries 
in the tree-dominated decays B — J/ipKs and in the penguin-dominated decays B — rj'Ks, which demands a deeper 
theoretical understanding. Besides, the branching ratios of the B — )■ r]'K decays predicted in the PQCD approach 
up to NLO are still much lower than the data. The value of f^, obtained in Sec. Ill A suggests a quantitative 
reexamination of the tree contribution from B — >■ rjcK to the direct and mixing-induced CP asymmetries and the 
branching ratios of the B^ — > ■q'-'^Ks decays. The Xcp factor for this study is defined as 

with the eigenvalue rjf — —1 and the weak phase /3. The direct and mixing-induced CP asymmetries are then derived 
from 
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For the i?° v'^s decays without the i]c mixing, one has 

where Vs are the CKM matrix elements, and T and P represent the tree and penguin amphtudes, respectively. 
These decay amplitudes have been evaluated up to NLO in the PQCD approach [34], which lead to Afip — 0.024 and 
^cW ~ 0.667^. The tree amplitude T,f^K has been also calculated in NLO PQCD 35], which gives the branching 
ratio Br{B^ — >■ Vr.K'^) = 5.5 x 10~^. Because the B meson decay constants fs — 0.21 GeV and fs = 0.19 GeV were 



adopted in [3j] and [35|, respectively, and the rjc meson decay constants /^^ = 0.478 GeV were derived in Eq. ([M)) . 
and /^^ = 0.42 GeV was adopted in [H, we multiply T^^k in [H by a factor (0.21/0.19)(0.478/0.42) = 1.26 for 
consistency. The increased NLO PQCD prediction Br{B'^ — )■ rjcK'^) = 8.7 x 10""* then agrees well with the data 
(8.7 ± 1.9) X 10~* [El]. Note that the relative strong phases among the above amplitudes T,fK, PifK and Tjj^k are 
known in the PQCD approach, so we do not encounter the difficulty mentioned in 22], and can derive the modified 
CP asymmetries without ambiguity. The inclusion of the B — > rjcK channel into Eq. (pg]) . 

yields A'^p = 0.023 and A™" = 0.664. It is seen that the rjc mixing causes a negligible effect on the CP asymmetries 
with A™]? moving slightly toward the central value of the data, 0.59 ±0.07 [5^. The result of A'j^p is consistent with 
the dat a 0. 05 ib 0.05 |52l |. Nevertheless, the rjc mixing brings the branching ratio Br{B° rj'K°) = 50 x 10"^ in NLO 
PQCD (ssj to 59 X 10~^, which becomes closer to the data (66.1 ± 3.1) x 10~^ [5^. Note that the enhancement of 
the above branching ratio due to the charm content of the 77' meson is larger than few percents estimated in [66j . 

Similarly, we investigate the impact of the rjc mixing on the B^ — > rjKs decays. Without the charm content of the 77 
meson, the NLO PQCD analysis gave = -0.128 and A'g'jS = 0.659 |i]- The rjc mixing then modifies the above 
values into A^'^ = -0.123 and A™]? = 0.644, namely, with a negligible effect. The branching ratio Br{B^ — > rjK^), 
becoming 2.3 x 10~^ from 2.1 x 10^^, is almost not changed by the rjc mixing. The result is a bit higher than the 
data Br{B° -J> rjK'^) = (1.12tQ 28) ^ lO"*" However, if using = -11° in the present work, the destructive 

interference for the B — > rjK decays from the rj-rj' mixing would be stronger, which will lower their branching ratios. 



IV. SUMMARY 



In this paper we have extended the rj-rj'-G mixing formalism constructed in our previous work to accommodate 
the rjc meson in a tetramixing scheme. The additional mixing angle between G and rjc was determined to be about 
11° from the observed widths of the rjc meson decays into light hadrons and 77. This mixing would have a leading 
impact on the rjc — !• gg width instead of on the rjc — > 77 one, such that the 0{asV^) corrections in NRQCD [i^ can 
be parametrized out. More precise measurement of the 77c total decay width and its decays into 77 can provide better 
constraints on the mixing scheme. Our tetramixing matrix was found to be consistent with that constructed from the 
SO(3) (X" 5*0(3) parametrization with a fit to data of relevant transition form factors. Contrary to general opinions in 
the literature, the present work suggests an reexamination of effects from the charm content of the 77''^ meson on the 
B — >■ rj'^'^K decays. 

We have shown that such a tetramixing scheme does increase the pseudoscalar density nigq to above the pion mass, 
which thus enhances theoretical predictions for the B —>■ rj' transition form factors and the B — )■ rj'K branching 
ratios [28]. It has been observed that the charm content of the 77' meson provides 18% enhancement of the B — rj'K 
branching ratios. The combined mechanisms can push the predicted values in NLO PQCD to the data easily, so the 
puzzle due to the large Br{B rj' K) is resolved. With this work, we postulate that the charm content of the rj' 
meson plays a more important role than the gluonic content in accommodating the large Br{B — >■ rj'K). Nevertheless, 
the 77c mixing has negligible effects on the direct and mixing-induced CP asymmetries of the B — > rj^'^^Kg decays, 
though the mixing-induced CP asymmetries move slightly toward the central value of the data. 

We have also investigated the impact of the tetramixing on the present theoretical and experimental observations 
of charmonium magnetic dipole transitions, and similar improvement is also seen: the gluonic content of the rjc 
meson decreases the decay widths T{J/ip,ip' — > 77/c) predicted by the nonrelativistic potential models by 7%. It has 
been confirmed that the rjc mixing does not modify the prediction in [25| for the pseudoscalar glueball mass in the 



^ The quoted value of A^^p differs from that presented in |34 ] . since the input of the weak phase f) has been corrected. 
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vicinity of 1.4-1.5 GeV. This result makes the 77(1405) meson an interesting candidate for the pseudoscalar gluebalL 
However, one should be aware of the complexity of underlying dynamics, such as the octet-glueball coupling |51j and 
intermediate meson rescattering in this mass region [63| , which certainly affect the determination of the pseudoscalar 
glueball mass. 
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